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The magnetic and thermodynamic properties of the two-dimensional (2D) quantum Heisenberg 
antiferromagnet (QHAF) that incorporates both a Dzyaloshinskii-Moriya (DM) and pseudo-dipolar 
(XY) interactions are studied within the framework of a generalized nonlinear sigma model (NLSM). 
We calculate the static uniform susceptibility and sublattice magnetization as a function of tem- 
perature and we show that: i) the magnetic-response is anisotropic and differs qualitatively from 
the expected behavior of a conventional easy-axis QHAF; ii) the Neel second-order phase transition 
becomes a crossover, for a magnetic field B _L CuCh layers. We provide a simple and clear expla- 
nation for all the recently reported unusual magnetic anisotropies in the low-field susceptibility of 
La2Cu04, L. N. Lavrov et at, Phys. Rev. Lett. 87, 017007 (2001), and we demonstrate explicitly 
why La2Cu04 can not be classified as an ordinary easy-axis antiferromagnet. 

PACS numbers: 74.25.Ha, 75.10.Jm, 75.30. Cr 



I. INTRODUCTION 

The usual starting point for the description of the mag- 
netism in undoped cuprates is the 2D isotropic Heisen- 
berg model 



H — J ^ Sj ■ Sj, 



(1) 



with only the AF super-exchange J between the nearest- 
neighboring Cu ++ spins, Sii Remarkably, a long- 
wavelength NLSM description of such model has been 
shown to correctly capture many of the features ob- 
served experimentally in the high-temperature param- 
agnetic phasei 2 ^ At lower temperatures, however, small 
anisotropies like the DM and XY interactions become im- 
portant, and the commonly used isotropic description is 
no longer adequate to describe the long-wavelength mag- 
netic properties of La2Cu04. In a recent paper, Lavrov 
et al. reported on some unusual anisotropies in the 
low-field magnetic susceptibility observed in detwinned 
La2CuC>4 single crystals^ Among the unexpected fea- 
tures one could mention: i) an almost featureless trans- 
verse in-plane susceptibility, Xa, for an extended temper- 
ature range; ii) the increase of the longitudinal suscepti- 
bility, Xb, as one approaches the Neel transition tempera- 
ture, Tjv, from the ordered side; iii) the anomalous T = 
hierarchy Xc > Xb > Xa, where a,b,c represent the low- 
temperature orthorhombic (LTO) crystallographic direc- 
tions of Fig. ^ 

In this article we investigate theoretically the magnetic 
response and thermodynamic properties of the square- 
lattice QHAF that incorporates both the DM and XY 
anisotropy terms by investigating the generalized NLSM 
corresponding to this microscopic problem. As it has 
been stressed in the context of one-dimensional systems^ 



the DM term generates an effective staggered magnetic 
field proportional to the applied uniform field. In the 
two-dimensional case the continuum field theory allows 
one to reproduce straightforwardly the same effect jS*l. 
In Refii the attention was focused on the (classical) spin 
configurations expected in an external field. In this paper 
we investigate instead the thermodynamical properties of 
the system, which enter in a crucial way in the evalua- 
tion of the static uniform susceptibility, x> as a function 
of temperature, T. In particular, we show that the cou- 
pling induced by the DM term between the magnetic field 
and the antiferromagnetic order parameter is particularly 
relevant when the order-parameter fluctuations become 
critical, leading to an anomalous magnetic response of the 
system, as observed for example hA As an outcome of 
our studies we observe that, while the isotropic nonlinear 
sigma model (NLSM) is well established to give a proper 
description of the magnetism in the high-temperature 
paramagnetic phase, T ^> for the ordered (bro- 

ken) phase, T < Tjy, we find more appropriate the adop- 
tion of a soft version of such fixed-length constraint, in 
order to correctly capture the physics of the longitudi- 
nal fluctuations. For the sake of clarity, we refer explic- 
itly to the La2Cu04 system throughout the paper, but 
the results presented here apply to any 2D square-lattice 
QHAF where the DM and/or XY anisotropy interactions 
are present, as for example other cuprates and nickelates 
La 2 Ni0 4 £ 
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II. THE LONG-WAVELENGTH LIMIT 



and 



We consider the following square-lattice single-layer 
S = 1/2 Hamiltonian for the La2Cu04 system 



H 



(id) 



{Si 



where Dy and T ^ are, respectively, 



Si ■ r ij ■ Sj , 

<M> (i,3) 

(2) 

the DM and 

XY anisotropic interaction terms that arise due to 
the spin-orbit coupling and direct-exchange in the low- 
temperature orthorhombic (LTO) phase of La2Cu04^ 
The direction and the alternating pattern of the DM vec- 
tors, shown in Fig. are a direct consequence of the 
tilting structure of the oxygen octahedra and of the sym- 
metry of the La2CuC>4 crystal. Throughout this work 
we adopt the LTO (bac) coordinate system of Fig. for 
both the spin and lattice degrees of freedom, and we use 
units where H = kn = 1. Thus we have that 



1 



= -j={-d,d,0), 



T> lk = -^=(d,d,0), (3) 
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where ij and ik label the Cu ++ sites (see Fig.^l, and d 
and T12.3 > are of order 10 -2 and 10 -4 , respectively, 
in units of j£ 

To construct the long wavelength effective theory for 
the above Hamiltonian we follow the standard procedure: 
we write 



Si(r) 
S 



^■ X 'n(x l ,r)+L(x l ,r), 



(4) 



where n and L are, respectively, the staggered and uni- 
form components of the spin and Q = (ir, tt). We then 
integrate out L and we obtain a modified NLSM where 
additional terms appear due to the DM and XY interac- 
tions (/? = 1/T) and 



S = \ dr / d 2 x { (d T n) 2 + c 2 (Vn) 2 + (D+ • n) 2 



(5) 



and the fixed length constraint n 2 = 1 is implicit (see also 
Refifiii). Here go is the bare coupling constant, related to 
the spin- wave velocity, cq, and renormalized stiffness, p s , 
through p s = co(l/ATgo — A/47r) (Si where A is a cutoff for 
momentum integrals (we set the lattice spacing a = 1) 
and N — 3 is the number of spin components. Finally, 



D+ = 2S(T>ij +Difc), 



see Fig. and 



r c = 32js 2 (r! - r 3 ) > o, 



(6) 



(7) 



and we neglected a small anisotropic correction to c com- 
ing from the XY term, since <C J. 

Observe that in the NLSM © only the sum of the DM 
and XY terms on the (ij) and (ik) bonds appears. This 
has several consequences: (i) the XY term generates a 
contribution to (jjj proportional to — r Qa n 2 , a — a,b,c, 
where T a p = (l/2)(Ty +T ik ) af) = diag(ri, Ti, r 3 ). Since 
n 2 = 1, such contribution reduces to a constant shift in 
the classical energy plus the T c n 2 term in Eq. (§J. This 
allows us to identify the ab as the first easy plane. The 
property that only the combination Ti — T3 affects the 
physics of the model (2J was an outcome of the numer- 
ical results of the RPA improved spin-wave analysis of 
Ref£, but no explanation for this effect was provided, 
(ii) The DM distortion enters the action (0 only through 
the vector D + = -D+x Q (x a ,X(,, and x c are the LTO 
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FIG. 1: Left: the hatched circles represent the 0~~ ions 
tilted above the &1O2 plane; the empty ones are tilted below 
it; small black circles are Cu ++ ions; bac orthorhombic coordi- 
nate system. Right: Schematic arrangement of the staggered 
magnetization (small black arrows) and DM vectors (open 
arrows). Center: continuum definition of D+. 



unit vectors), rendering be a second easy plane. When 
the system orders we find the staggered magnetization 
at (n) = (ToXb, as observed experimentally^ because the 
orientation along x a or x c would cost an energy D\ or T c , 
respectively. Moreover, since the uniform magnetization 
is given by^i 



(L> = — «n) xD+), 



(8) 



we find that, in the AF phase, (L), is directed along x c , 



3 



so that the Cu ++ spins are canted out of the ab plane 
and confined to the be plane. 



III. MAGNETIC SUSCEPTIBILITY 



To evaluate the static uniform susceptibility, Xa, we 
rederive from the microscopic Hamiltonian the quan- 
tum action (0 in the presence of a magnetic field B (in 
units of gsHB/fr — 1, where gs ~ 2 is the gyromagnetic 
ratio and /is is the Bohr magneton)^! 



5(B) = S(d T n -> d T n + iB x n) 



— / dr / d 2 x B ■ (D + x n). 

go co Jo J 



(9) 



The last term in the above equation, which couples the 
DM vector D + and the staggered magnetization n, is 
responsible for many of the unusual features observed 
experimentally in La2Cu04. However, this coupling is 
generic to any square-lattice system, as other cuprates 
and nickelates, where symmetry guarantees that the DM 
vectors alternate in sign between neighboring bonds £ 

The zero-field magnetic susceptibility is calculated in 
linear response theory as 



1 d 2 \oeZ 



0V dBl 



(10) 



B=0 



where Z(B) — J 2?ncxp {— 5(B)} is the Euclidean par- 
tition function for the action ©, and we obtain 



Xa 







Xb = Xt + 



Xc 



go co 
go co 



-Gc(O.O), 



where 



G Q 1 (k, LO n 



goco 5oCo 



= c^k 2 



G b (0,0), 



Mi 



(11) 



(12) 



is the inverse propagator for the magnetic modes a, 6, c, 
with gaps M a and 



Xc 



1 

w 



J2 { Gp (q) + G 7 (q) - AloI Gp (q) G 7 (q) } , 



q=(k,, 



(13) 

is the traditional "uniform" contribution to the suscepti- 
bility, with (a/37) = {obc) and its permutations. 

The set of Eqs. ijTTll . following from the quantum ac- 
tions © and ©, is the main results of this article. A 
few remarks are in order now concerning the usefulness 
of the NLSM approach to the model <|2j with respect to 
other approaches 1 W' 11 First, the action © allows one to 
easily track down the source of the anisotropic magnetic 
response, as coming from the term B • (D+ x n). In fact, 



depending on the directions of D + and of the applied 
field, there will be (or not) additional terms proportional 
to G a in Eq. itTTll . responsible for the unusual magnetic 
response. Second, this coupling is generic to any system 
where the lack of inversion-center symmetry allows for 
an oscillating DM interaction between neighboring spins. 
Third, Eqs. l(TT|l are exact, within linear response theory, 
even though the evaluation of the zero-field thermody- 
namic quantities <Jq{T) and M a (T) will require, in gen- 
eral, a certain degree of approximation. Nevertheless, we 
show now that all the qualitative features observed ini 
are already present at the mean-field level. 

The thermodynamic properties of the model (2J can 
be investigated by means of a large-iV expansion^ At 
N = 00, and for T > Tjy, we find ctq = and 



Ml 



4c 



M 



r2 



Dl- 



M: 



r c+c 2 r 2 , (14) 



where the correlation length, £, is calculated through an 
averaged constraint equation 1 = NI±(£). Here I± = 
(l/2)(J a + / c ) with 

T goco sr^ n n \ SoT ( sinh (c A/2T) 1 
" = Wj^ GaiK ^ = ^0 bg ( sinh (MJ2T) ) • 

(15) 

Since both transverse spin-wave modes, a and c, are 
gapped, the 2D system orders at a finite Neel temper- 
ature, Tn, defined by 1 = NI±(£ = 00)^ 

For the ordered phase, T < Tjv, instead of replacing 
the local constraint by an average one, we introduce a 
potential of the type (wo/2)(n 2 — l) 2 in the action (§J, 
making the constraint locally soft. In this case we obtain 
a generalized linear sigma model (LSM), which leads to 
an ordered phase described by: 



1-NI±, Mi 



2cr 2 U , 



Mi = D±, Mi 



(16) 

We should emphasize that the only new parameter here 
is the energy scale uo — (7J) 2 , where 7 is a fitting pa- 
rameter that will be fixed later through comparison with 
experiments. However, an inspection of Eq. 11 (ill , in- 
dicates that neither Tjv nor the values of <tq and M a , c 
depend on uq. The only difference is the physics of the 
longitudinal fluctuations, for T < Tjy. While the hard 
constraint yields a longitudinal mode always gapless, 13 
the soft one gives a mass for the longitudinal mode be- 
low T/v that is proportional to the strength of the order 
parameter. 

We fix the large- N parameters for La2Cu04 by com- 
paring the correlation length, £, with the experimental 
data in the paramagnetic phase, and we use M a (T = 
Q)=D+ = 2.5 meV and M C (T = 0) = ^JT~ C = 5.0 meV, 
for the zone center in-plane and out-of-plane spin-wave 
gaps as determined from neutron scattering^ In the in- 
set (a) of Fig.[2]we compare £ -1 obtained for p s = 0.1J, 
c = 1.3J and J — 100 meV with the curve Cxp, which 
represents the best fit of the experimental data at T > T, 



N 



obtained in the isotropic NLSMl*i (so that £ e x P deviates 
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FIG. 2: (Color online) Temperature dependence of the mag- 
netic susceptibilities Xa> Xb and Xc from Hill . Insets: (a) 
temperature dependence of M a ,b,c and (b) 



from the experimental data near Tjv). When l/N cor- 
rections to our calculations are included one expects that 
the ratio Tn/J is reduced^ which allows for larger and 
more realistic values for p s , Co, and Nevertheless, 
with the above parameters we find ao(T = 0) w 0.46 
which leads to an effective moment of the Cu ++ spins, 
M = SsSagfiB ~ 0.46/is, in good agreement with the 
experiments^ The resulting M a> b c as a function of tem- 
perature are shown in the inset (b) of Fig.|2]for 7 = 0.15, 
a choice that will be justified below. 



tween the decrease of \a and the increase of the term 
We turn now to Xc and Xb- As T — > Tjv, 



1 



ffoco goco M b z 



(17) 



diverges because of the vanishing of the mass of the lon- 
gitudinal mode, M b (T — > TjJ) — > 0, and this is associated 
with the ferromagnetic ordering of the canted momentsiii 
It is essential to have a finite Mb below T/v in order to cor- 
rectly reproduce the behavior of Xc in the ordered phase. 
So, although the hard-spin model correctly reproduces 
the vanishing of M& as T — > T^, it is not adequate for 
the ordered phase, where Mb — leading to a diverging 
Xc for all T < Tn ■ That is why we have adopted the soft 
version of the constraint below Tn- 

Finally, we observe in Fig.|2lthat Xb is not divergent at 
Tjv but exhibits a well pronounced peak. This is a conse- 
quence of the term Z?+G c (0, 0)/goc o = D^_/(g c M 2 ). 
According to Eq. ll"14t . as we approach Tn from the 
paramagnetic side t]~ 2 goes to zero and D^/(gocoM 2 ) 
approaches the value D^_/(g coT c ) « l/(4g co), which is 
controlled by the ratio between the DM and XY interac- 
tion terms. Observe that the absence of a similar feature 
in Xa is due to the fact that, for B || a, B- (D + x n) = 0, 
so Xa has only the uniform contribution Xa- Thus, the 
increase in Xb as T — > T^ is clearly a result of the sum of 
its uniform part, Xb > which is monotonically increasing 
with T, and the term D+/(goCoM 2 ), which is constant 
below Tjv and decreases instead as ~ £ 2 above Tjv. 



A. Temperature dependence of Xa,Xb and Xc 

Let us now discuss how all the unusual features of 
the susceptibilities reported in 4 follow straightforwardly 
from Eqs. l(TT|l already at the mean-field (large N) 
level (see Fig. EJ. At large temperatures, T ^> Tn, 
all three Xa.b.c exhibit the usual linear-T behavior ob- 
served experimentally^ and expected for the uniform 
susceptibility of the isotropic NLSMi^ As T — > 0, on 
the other hand, Xb while both x" and Xa saturate 
at very small, but nonzero, values. The smooth evolu- 
tion through the transition experimentally observed in 
Xa is due to the fact that Xa contains terms like I a in 
Eq. II 1511 , which are regular for the transverse modes (al- 
ways gapped by the DM and XY interactions) and di- 
verge logarithmically at Tn for the longitudinal mode 
as Mb(T — > Tjv) —> 0. This divergence appears as a 
hump in the numerical evaluation of Xa presented in Fig. 
El and is an artifact of both the low dimensionality and 
of the large- N limit (a coupling between CuC>2 layers or 
a nonzero anomalous dimension r\ for the propagator of 
the longitudinal mode regulates such logarithmic infrared 
divergence of lb at T N and wipes out the feature in Xa)- 
Finally, the weak increase observed in Xa as T — > 0, also 
observed experimentally, is due to the compensation be- 



B. Unusual x«(0) < Xb(0) < Xc(0) hierarchy 

A comment is in order now concerning the hierarchy of 
the T = susceptibilities in Fig. [21 For an ordinary easy- 
axis AF with different in-plane and out-of-plane trans- 
verse gaps, M a < M C) it is expected that the T = 
uniform susceptibilities should satisfy Xbify < X"(0) < 
X"(0), with the lowest value for the easy axis and with 
a transverse susceptibility smaller in the direction of the 
smaller gap. This is in fact the result that one obtains 
after dropping the terms proportional to D\ in Eq. itTTll . 
However, the unusual coupling between B and n in © 
gives rise to the terms proportional to D\ in x&(0) and 
Xc(0), which, in turn, lead to the following T = values 
for the three susceptibilities 

~_£o_ -J_^± ~ 1 D\ 

go Co 5o Co M£ 50 co M b z 

(18) 

resulting in the unusual hierarchy of the zero- 
temperature values and rendering La2Cu04 an example 
of an unconventional easy-axis antiferromagnet. Observe 
that Xc(0)/Xa(0) > 1 always, and we use the T — in- 
tercept of Xc to fix 7 = 0.15. Finally, for our choice 
of parameters the ratio Xb(0)/xa(0) ~ 1.1, in relatively 
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good agreement with the experiments^ but it depends in 
general on the values of the transverse massesiAi 







Exp x 




X c — 




/ Ac 








0.2 

T/J 



FIG. 3: (Color online) Plot of a , M b , and M^ p from Eqs. 
12111 and 12011 . M^ xp is fitted to the experimental data for the 
inverse correlation length above Tat, see discussion in the text 
below Eq. 12111 . Inset: Xc in units of emu/g and experimental 
points from 4 . \T V is obtained from Eqs. Ill IP using M£ xp . 
Here p s = 0.07J, 7 = 0.5, a — 5.3 A, the out-of-plane lattice 
parameter d — 13.5 A, and we added a T-independent shift 
of x v v ~ x 10~ 7 emu/g to account for the van Vleck 
paramagnetic susceptibility!™ 



IV. EFFECS OF A FINITE MAGNETIC FIELD 
B 1 C11O2 LAYERS 

The result for \c plotted in Fig. |2can be further im- 
proved by considering explicitly the nontrivial effects of 
the combination B _L ab plane and D + ^ on the ther- 
modynamic properties of the theory ©. In fact, in this 
case the last term in Eq. Q becomes — (1/goCo) / h rib, 
where h = |B X D+|, Within the NLSM formulation the 
averaged constraint equation defines a self-consistency 
equation for the order parameter given by: 



1- NI ± (Ml + h/a ) 



(19) 



where the masses appearing in the transverse spin fluc- 
tuations 1151 are themselves a function of oq : 



Ml = Dl 



h 



Mi 



r, 



h 

Co 



B 



(20) 



As one can see, besieds a (quantitatively small) hard- 
ening of the c mode due to the applied field, we find 
that the term h/ao plays a a role similar to the corre- 
lation length, c 2 £~ 2 , in Eq. 11411 . The result is that <jq 
never vanishes for h ^ (see Fig. |2J , and because of the 
DM interaction, a nonvanishing B _L ab plane transforms 
the Neel second order phase transition into a crossover. 
In analogy with the discussion of the case B = 0, the 



proper description of the longitudinal-field fluctuations 
below the crossover temperature can be done within the 
LSM, where the order-parameter equation and the longi- 
tudinal mass read: 



u a (al-l+NI ± {Ml+h/a ,T)) = h, M 2 = 2u a 2 + — . 

CTO 

(21) 

Observe that the crossover temperature obtained from 
Eq. iflfll) and Eq. l(2"l")l is almost the same. To inter- 
polate between the two solutions we define (M^ xp ) 2 — 
max(h/ <Jo, 2uo<Jq + h/ao) and we fix the parameter val- 
ues by fitting M^ xp to the the experimental data on the 
correlation length above Tjy. As far as the c-axis suscep- 
tibility is concerned it is now evident from the inset of 
Fig-GDthat no i on g er diverges (Mb is always finite) but 
becomes peaked at a crossover temperature at which Mb 
has a minimum (see Fig.EJ. Observe also that for \a and 
Xc the effects of a finite magnetic field are instead negli- 
gible, leading only to small quantitative differences with 
respect to the calculation presented before for B = 0. 



V. CONCLUSIONS 

In conclusion, we have derived the uniform magnetic 
susceptibility within the long- wavelength effective theory 
©-© for the single-layer square-lattice QHAF with DM 
and XY interactions l|2j . We obtained that the magnetic 
response is anisotropic, in remarkable agreement with the 
experiments and differs from the expected behavior 
for a more conventional easy-axis QHAF. Due to the pres- 
ence of the DM term, the uniform magnetic field gener- 
ates an effective staggered field proportional to both the 
DM interaction and to the applied field. We showed that 
the coupling between the magnetic field and the stag- 
gered order parameter is particularly relevant when the 
order-parameter fluctuations become critical, leading to 
an anomalous magnetic response of the system, as ob- 
served for example in-. Moreover, this same coupling 
is responsible for the unusual zero temperature hierachy 
of the susceptibilities, rendering La2Cu04 an example 
of an unconventional easy-axis antiferromagnet. Finally, 
we considered the effect of a finite magnetic field on the 
Xc susceptibility. We found that, for a nonzero B _L ab 
plane, the Neel second order phase transition becomes a 
crossover, see Fig. in analogy with the beahvior of a 
ferromagnet in a finite magnetic field. The finite-field ef- 
fects are instead irrelevant as far as the susceptibilities in 
the a and b direction are considered. However, the case of 
B || b presents interesting outcomes as far as the selection 
rules for one-magnon Raman scattering are concerned, 
as it has been discussed recently in Ref^. While Eq. 
and mill have been derived for the case of La2Cu04 
materials, the present analysis can be extended to other 
QHAF systems where the lack of inversion-center sym- 
metry allows for an oscillating DM interaction between 
neighboring spins. 



6 



VI. ACKNOWLEDGEMENTS tro Neto, and R. Gooding. We also thank A. Lavrov, for 

fruitfull discussions and for providing us with the exper- 
The authors have benefitted from discussions with imental data for X c from*, shown in Fig. (inset). 
A. Aharony, Y. Ando, B. Binz, C. Castellani, A. H. Cas- 



* Electronic address: barbosa@phys.uu.nl 

1 M. A. Kastner, R. J. Birgeneau, G. Shirane, and Y. Endoh, 
Rev. Mod. Phys. 70, 897 (1998). 

2 S. Chakravarty, B. I. Halperin, and D. R. Nelson, Phys. 
Rev. B 39, 2344 (1989). 

3 A. V. Chubukov, S. Sachdev, and J. Ye, Phys. Rev. B 49, 
11919 (1994). 

4 A. N. Lavrov, Y. Ando, S. Komiya, and I. Tsukada, Phys. 
Rev. Lett. 87, 017007 (2001). 

5 M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79, 2883 
(1997); Phys. Rev. B 60, 1038 (1999). 

6 L. Benfatto, M. B. Silva Neto, V. Juricic, and 
C. Morais Smith, to appear in the Proceedings of the Inter- 
national Conference on Strongly Correlated Electron Sys- 
tems - SCES05, Vienna, 2005. 

7 J. Chovan and N. Papanicolaou, Eur. Phys. J. B 17, 581 
(2000). 

8 L. Shekhtman, O. E. Wohlman, and A. Aharony, Phys. 
Rev. Lett. 69, 836 (1992); W. Koshibae, Y. Ohta, and 
S. Maekawa, Phys. Rev. B 50, 3767 (1994). 

9 K. V. Tabunshchyk, and R. J. Gooding, Phys. Rev. B 71, 
214418 (2005). 

10 D. Vaknin, S. K. Sinha, D. E. Moncton, D. C. Johnston, 
J. M. Newsam, C. R. Safinya, and H. E. King Jr., Phys. 
Rev. Lett. 58, 2802 (1987). 

11 T. Thio, T. R. Thurston, N. W. Preyer, P. J. Picone, 



M. A. Kastner, H. P. Jenssen, D. R. Gabbe, C. Y. Chen 
and R. J. Birgeneau, and A. Aharony, Phys. Rev. B 38, 
R905 (1988); T. Thio, and A. Aharony, Phys. Rev. Lett. 
73, 894 (1994). 

12 The interlayer coupling, J±, plays also an important role 
for the Neel transition in La2Cu04. We will nevertheless 
proceed with our 2D model since it already captures the 
essential physics reported inA 

13 V. Yu. Irkhin, and A. A. Katanin, Phys. Rev. B 57, 379 
(1998). 

14 B. Keimer, R. J. Birgeneau, A. Cassanho, Y. Endoh, 
M. Greven, M. A. Kastner, and G. Shirane, Z. Phys. B 
91, 373 (1993). 

15 B. Keimer, N. Belk, R. J. Birgeneau, A. Cassanho, 
C. Y. Chen, M. Greven, M. A. Kastner, A. Aharony, Y. En- 
doh, R. W. Erwin, and G. Shirane, Phys. Rev. B 46, 14034 
(1992). 

16 T. Nakano, M. Oda, C. Manabe, N. Momono, Y. Miura, 
and M. Ido, Phys. Rev. B 49, 16000 (1994). 

17 In the RPA improved spin- wave analysis o& Xa(0) ~ 
l/poco and as a consequence Xb(0) > Xa(Q) is only obtained 
if M a > M c , in contrast to the experimental situation 1 *. 

18 M. B. Silva Neto and L. Benfatto, |cond-m at/0507103 to 
apper in Phys. Rev. B (2005). 



